Entanglement and decoherence: fragile and robust entanglement 



O 
(N 



> 
oo 

m 
rn 

rn 
O 



X 



Jaroslav Novotny,^'^ Gemot Alber,^ and Igor Jex^ 

^Department of Physics, FNSPE, Czech Technical University in Prague, 
Bfehovd 7, 115 19 Praha 1 - Stare Mesto, Czech Republic 
^ Institut fiir Angewandte Physik, Technische Universitdt Darmstadt, D-64289 Darmstadt, Cermany 

The destruction of entanglement of open quantum systems by decoherence is investigated in the 
asymptotic long-time limit. Starting from a general and analytically solvable decoherence model 
which does not involve any weak-coupling or Markovian assumption it is shown that two funda- 
mentally different classes of entangled states can be distinguished. Quantum states of the first class 
are fragile against decoherence so that they can be disentangled asymptotically even if coherences 
between pointer states are still present. Quantum states of the second type are robust against de- 
coherence. Asymptotically they can be disentangled only if also decoherence is perfect. A simple 
criterion for identifying these two classes on the basis of two-qubit entanglement is presented. 
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Every natural physical quantum object is in contact 
with an environment. A consequence of the resulting in- 
evitable interaction of such an open quantum system with 
its environment is decoherence [l|-|3| . For sufficiently long 
interaction times this process leads to an environment- 
induced selection of so-called pointer states which re- 
main at least approximately undisturbed while any of 
their superposition states decays quickly into a mixture. 
This asymptotic loss of quantum coherence is responsible 
for the appearance of classical features of open quantum 
systems [l|-|3[. Apart from that decoherence is also a 
process which tends to destroy entanglement. Neverthe- 
less, the recently discussed time-dependent phenomenon 
of sudden death of entanglement [4| has already hinted 
at some subtle relations between decoherence and disen- 
tanglement in the time domain. Thereby, it has been 
observed that entanglement of an open quantum system 
may already disappear after a finite interaction time with 
its environment while complete decoherence requires an 
infinite amount of time. A further exploration of the 
subtleties between entanglement and decoherence is not 
only of fundamental physical interest but is also of sig- 
nificance for practical applications in the area of quan- 
tum information science |5| in which entanglement is a 
key resource. For the realization of large quantum infor- 
mation processors, for example, it is of vital importance 
to understand how decohcrence-induced loss of entangle- 
ment of an open quantum system scales with the size of 
an environment in the limit of arbitrarily long interac- 
tion times. Under which conditions are already finitely 
sized environments capable of disentangling open quan- 
tum systems completely? Are there classes of entangled 
quantum states whose properties differ significantly un- 
der decoherence? It is a main purpose of this letter to 
explore these questions. 

For this purpose a sufficiently general model of deco- 
herence is needed which is exactly solvable without fur- 
ther simplifying assumptions concerning, e.g., sizes of en- 
vironments, initially prepared quantum states, interac- 



tion strengths, or correlation times between system and 
environment. Of particular interest are generic asymp- 
totic long-time properties of such a general decoherence 
model which exhibit clearly the intricate interplay be- 
tween decoherence and destruction of entanglement. 

Motivated by the practical significance of elementary 
distinguishable two-level systems (qubits) for purposes 
of quantum information processing in the following we 
present a qubit-based class of such decoherence models. 
Thereby, we start from the recent observation that de- 
coherence of a single system qubit can be modeled by 
a sequence of 'controlled-U' unitary transformations in 
which the single system qubit acts as a control and the 
environmental qubits act as targets @. This result was 
obtained under the simplifying assumptions that the en- 
vironment is formed by an infinite number of qubits, that 
each qubit is prepared initially in the same state so that 
all environmental qubits are uncorrelated, that the envi- 
ronmental qubits do not interact among themselves, and 
that each environmental target qubit interacts with the 
single system qubit only once. In order to overcome all 
these restrictive approximations let us consider as a nat- 
ural generalization k system qubits which interact with n 
environmental qubits by a sequence of elementary inter- 
actions or 'collisions'. All of these 'collisions' are assumed 
to be well separated in time. In each of them a pair of 
qubits, i.e. a control qubit i and a target qubit j, is se- 
lected randomly and a controlled unitary transformation 
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is applied. Thereby, Ij is the unit operator acting on 
qubit j and 

uf^ = cos(A(|0),,(0| - |1),,(1|) -f sin0(|O),,(l| + |1),,(0|) 

(2) 

denotes the unitary one-qubit transformation acting on 
the target qubit j. Such controlled unitary couplings 
have already been investigated in the context of vari- 
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ous decoherence models [l|, |3|] because they do not af- 
fect the computational basis states |0)i and of the 
control qubit i and at the same time they can decrease 
any quantum coherence between these two system states. 
Within our generalized decoherence model each system 
qubit z is a possible control qubit for an elementary inter- 
action U^^^ with any environmental target qubit j. Fur- 
thermore, all environmental qubits also interact among 
themselves. However, in order to guarantee a well de- 
fined pointer basis in our decoherence model we assume 
that the system qubits do not interact among themselves. 
The interaction pattern characterizing which qubits can 
be coupled by an elementary interaction U^^^ between 
randomly chosen qubits i and j can be encoded in a con- 
venient way in a weighted and directed interaction graph 
. Its set of vertices represents the n + k qubits of the 
total qubit system and the set of directed edges E rep- 
resents the possible interactions among the qubits. Each 
of these edges e = ij G E is weighted with a probability 
Pf, with which qubits i and j are coupled by the unitary 
two-qubit transformation Ue'^\ These probabilities are 
normalized to unity, i.e. X^eeEPe = 1- A simple example 
of such a qubit network is depicted in Fig. [1] 




FIG. 1: Schematic representation of a 
system-environment qubit network: k = 2 qubits of the 
open quantum system (grey region) interact with n = 2 
environmental qubits. The directed edges indicate all 

possible elementary two-qubit couplings. The tail 
(head) of an edge symbolizes the control (target) qubit. 

Thus, within this decoherence model the quantum 
state p{N) resulting from N completed elementary in- 
teractions is changed by the (A^-|- l)-th elementary inter- 
action to the quantum state 

p{N + 1) = ^p,;7i^)p(7V);7i^)t ^ V{p{N)). (3) 

eGE 

The map is a random unitary operation Q and de- 
scribes an elementary interaction or 'collision' between 
system and environment. General properties of quantum 
states p{N) resulting from repeated applications of such 
quantum maps have been studied recently @. In par- 
ticular, it has been shown that in the asymptotic limit 
of ^ 1 the quantum state p{N) becomes independent 
of the probability distribution {pe,e £ E} and that it 
is determined uniquely by a linear attractor space. This 
latter space is formed by the maximal set of all possible 



orthonormal solutions Xx.i of the eigenvalue equations 

U^'f'^Xx,rU^^'^^ = XXx,r for all e e E (4) 

with |A| = 1. The index i distinguishes different mutu- 
ally orthonormal solutions with eigenvalue A. If these 
attractors are known, in the asymptotic limit ^ 1 the 
quantum state resulting from an initially prepared state 
Pin is given by 

p{N) = V^{p^n)= A^Tr{/}„,Xl_jlA,. (5) 

|A| = l/j 

With the help of the methods developed in Ref. Q 
it can be shown that within our decoherence model for 
n > 2 the only possible eigenvalue of Eqs.Q is given by 
A = 1 and its eigenspace is spanned by the 4*^ + 3 x 2^^ + 1 
linear independent solutions 

\x){x\ ® /„, |0)(a;| (g) |O„)(0„|, |a;)(Ofc| ® |(/)„)(0„|, 
\x){y\ ® |0„)(0„|, |0fc)(0fe| ® |0„)(0„| (6) 

with the pure quantum states |0„) = |0)®", |0) = 
cos((/)/2)|0) -|-sin(^/2)|l), |(/)„) = |(/))®". Vectors |a;) and 
\y), with x,y E 2'^, are arbitrary elements of the compu- 
tational basis of k qubits. Note that the quantum state 
is not affected by the unitary transformation u'^ of 
Eq.©, i.e. u-^l./)) = 10). 

In order to investigate the decoherence process affect- 
ing the k system qubits in the asymptotic limit of large 
numbers of elementary interactions A'^ ^ 1 let us consider 
an initially prepared quantum state pin = p^'^' ® p'^-* 
which does not contain any correlations between sys- 
tem S and environment E. Using the attractor space 
of Eq. ([U and tracing out the environment in Eq. ([5]) the 
asymptotic quantum state of the open quantum system 
S is given by 

/)(^)=diag(p(^))-K (^n\p^''^\c^n)&x){y\ 

(7) 

with diag{p^^^) denoting the diagonal part of the system 
state p'^^^ with respect to the pointer basis {\x),x G 2'"'}. 
At this point it is worth emphasizing that in view of the 
general results of Ref.0 some of the conditions leading 
to Eq.© can be relaxed. Thus, it turns out, for exam- 
ple, that Eq.© is already valid if at least one pair of 
environmental qubits is interacting. 

According to Eq.Q the influence of the environment 
E on the asymptotic quantum state of the system S is de- 
scribed by the decoherence factor < r = ((^„ |0„) < 
1 which depends only on environmental properties, such 
as its size n and the initially prepared quantum state 
p^-^\ In particular, three different cases can be distin- 
guished. If all environmental qubits are initially prepared 
in the eigenstate of the interaction of Eq.([2]), i.e. 
r = 1, the system state p^^^ remains unchanged. If the 
environmental state p^^^ has zero overlap with quantum 
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state |(^„), i.e. r = 0, we observe perfect decoherence in 
the asymptotic limit ^ 1. In all other cases asymp- 
totic decoherence results in an only partial suppression of 
coherences between the pointer states of the open system 
S. In particular, if the initial state of the environment 
is a factorized state ^"^"^ for example, the decoherence 
factor is given by r = Thus, coherences be- 

tween pointer state of the system S decrease exponen- 
tially with increasing number n of environmental qubits. 
In this case perfect decoherence can only be obtained for 
an infinite environment. Furthermore, if environmental 
qubits share initial correlations even in the case of an en- 
vironment with infinitely many qubits perfect asymptotic 
decoherence cannot always be achieved. Let us consider 
an initially prepared pure environmental state of the form 
Penv = \Xn){Xn\ with \xn) = COS a\<pn) +siiia\un) and 
{(^nli^n) — 0, for example. The associated decoherence 
factor is given by r = cos^(a) and is thus independent of 
the number of environmental qubits n. 

Let us now investigate how this decoherence process 
destroys entanglement in the asymptotic limit of many 
elementary interactions ^ 1 between system and en- 
vironment. In order to avoid problems with appropri- 
ate measures of entanglement let us focus on an arbi- 
trary two-qubit subsystem of an open quantum system 
with k > 2 system qubits and n environmental qubits. 
Thus, the entanglement of this two-qubit subsystem can 
be quantified by its concurrence fioj . 

Consider first of all an initially prepared pure state 
Pi = of these two system qubits with \ipi) = 

a\00) + b\ll) and |ap -|- |5p = 1. For an arbitrary envi- 
ronmental state p*^^^ the concurrence C of the resulting 
asymptotic quantum state (O is given by 

C(M^) ^ 2|a6|(</.„|p(^)|0„) = (8) 

Thus, the resulting asymptotic entanglement is equal to 
the initial degree of entanglement multiplied by the de- 
coherence factor r. Therefore, in the asymptotic limit 
A'^ ^ 1 the disentanglement of the originally prepared 
entangled state p[^^ = |V'i)(V'i| is governed by the same 
dependencies concerning the size n and the initially pre- 
pared state of the environment p*^^) as the decoherence 
factor r. 

Let us now consider an uncorrclatcd state of sys- 
tem and environment with the two system qubits ini- 
tially prepared in the maximally entangled state |'02) = 
l/2(-|00) + |01) + |10) + For an arbitrary envi- 

ronmental state p^^^ the concurrence C of the resulting 
asymptotic quantum state ([7]) is given by 

C'ipi'J,) = I max {O, 3(0„|p(^) - l} . (9) 

Thus, the resulting asymptotic entanglement is drasti- 
cally different from the previous case. As soon as the 
decoherence factor decreases below the threshold value 



of 1/3 the degree of entanglement drops to zero in the 
asymptotic limit A^ ^ 1. The entanglement of the open 
quantum system is now very sensitive to the size n and to 
the state p*-^) of the environment. In particular, consider 
an n-qubit environment initially prepared in a factorized 
quantum state p'-^^ = with ((A„|/5(^)|(/i„) < 1. This 
implies the relation {4>n\p^^^\4'n) = Therefore, 
as soon as the number n of qubits in the environment 
exceeds the critical value of 

nse,^\- (10) 

the asymptotic entanglement of the system state van- 
ishes, ([x] denotes the largest integer less or equal to x.) 
This is in extreme contrast to the previous case studied 
in Eq.® in which entanglement decreases exponentially 
with the size n of the environment and never vanishes 
completely for any finite size of the environment. 

The numerical simulations depicted in Fig. [2] show the 
dependence of the concurrence C of the two qubits ini- 
tially prepared in the quantum state |'!/'2) on the number 
A^ of elementary interactions between system and envi- 
ronment. The asymptotic limit is already achieved after 
a few elementary interactions. In the exceptional case 
of an environment of minimal size, i.e. rt = 1, Eqs.Q 
also allow a solution with eigenvalue A ~ — 1 which re- 
sults in a non-stationary asymptotic limit leading to an 
oscillatory dependence. In all other cases the asymp- 
totic entanglement is stationary. In the case depicted in 
Fig. [2] the critical size of the environment is given by 
nsep = 4 so that for all smaller sizes of the environment 
asymptotically the two-qubit system is not disentangled 
completely. 




FIG. 2: Dependence of the concurrence C of the 
two-qubit state p^^^^ (N) on the number of elementary 
interactions N and on the size of the environment n 
(n = 1 full, n = 2 dotted, n = 3 dashed-dottcd, n = 4 
dashed): The initial state of the environment (system) 
is given by |,5)^" with \S) = sin(7r/6)|0) + cos(7r/6)|l) 

(1-02)) and u^"^) is given by Eq.© with </) = 27r/3. 

Thus, with respect to asymptotic decoherence- induced 
disentanglement of an open two-qubit system two classes 
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of initially prepared entangled states p^"^) can be distin- 
guished. The first class consists of quantum states which 
are fragile against decoherence in the sense that their 
asymptotic entanglement vanishes already for non-zero 
decoherence factors r (compare with Eq.®). In contrast, 
quantum states p^^^ of the second class are robust against 
decoherence-induced disentanglement in the sense that 
the two qubits can be disentangled only if decoherence 
is perfect, i.e. r = (compare with Eq.dH)). As a result 
these two classes of initially prepared two-qubit states ex- 
hibit a very different behavior with respect characteristic 
features of the decoherence-inducing environment. 

Within our general decoherence model even a simple 
necessary and sufficient condition can be derived whether 
the entanglement of an initially prepared two-qubit sys- 
tem state p'-'^^ of an open system is fragile or robust. In 
view of the general form of the asymptotic system state 
of Eq. d?!) the formal arguments presented in Ref. [lH lead 
to the conclusions that a two-qubit state p^^^ is fragile 
with respect to environmental decoherence if and only if 
its density matrix elements in the pointer basis satisfy 
the relation 



two-qubit subsystems of an open qubit system a simple 
criterion for robustness and fragility has been presented. 

The very existence of fragile and robustly entangled 
states in open quantum systems is a consequence of the 
characteristic structure of the asymptotic quantum states 
described by Eq.(9). Without the restriction < r < 1 
on the decoherence factor a similar structure also appears 
in other recently discussed decoherence models [l| 
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Therefore, it can be conjectured that the existence of 
fragile and robustly entangled quantum states is a general 
phenomenon accompanying any decoherence process. In 
view of their significantly weaker sensitivity to decoher- 
ence robustly entangled quantum states are expected to 
play a particularly significant role in the further develop- 
ment of quantum information processing and in its efforts 
to push entanglement as a characteristic quantum phe- 
nomenon as far as possible into the macroscopic domain. 
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(00|p(^)|00)(10|/5(^)|10)(01|/5(^)|01)(ll|p(^)|ll) ^ 0. 

(11) 

In summary, a class of decoherence models has been 
presented by which the intricate interplay between de- 
struction of entanglement and decoherence resulting from 
the asymptotic long-time interaction of an open qubit 
system with its environment can be investigated analyti- 
cally. It does not involve any further simplifying assump- 
tions, such as the ones involved in Markovian models or 
in models in which each environmental qubit interacts 
only once with a system qubit. Within this framework it 
has been shown that two characteristic classes of entan- 
gled quantum states can be distinguished, namely fragile 
and robustly entangled states. They exhibit significantly 
different dependencies on environmental properties. For 



[1] E. Joos et al., Decoherence and the Appearance of a Clas- 
sical Worldin Quantum Theory, (Springer, Berlin, 2010). 

[2] M. A. Schlosshauer, Rev. Mod. Phys. 76, 1267 (2005). 

[3] W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 

[4] T. Yu and J. H. Eberly, Phys. Rev. B 68,165322 (2003). 

[5] M. A. Nielsen and I. L. Chuang, Quantum Computation 
and Quantum Information (Cambridge UP, Cambridge, 
2000). 

[6] M. Ziman and V. Buzek, Phys. Rev. A 72, 022110 (2005). 

[7] J. Novotny, G. Alber, I. Jex, NJP, in print (2011). 

[8] I. Bengtsson and K. Zyczkowski, Geometry of Quantum 

States (Cambridge UP, Cambridge, 2006). 
[9] J. Novotny, G. Alber, I. Jex, Cent. Eur. J. Phys. 8, 1001 

(2010). 

[10] S. Hill, W. K. Wootters, Phys. Rev. Lett. 78, 5022 (1997). 
[11] J.-H. Huang, S.-Y. Zhu, Phys. Rev. A 76, 062322 (2007). 
[12] F. M. Cucchietti, J. P. Paz, W. H. Zurek, Phys. Rev. A 
72, 052113 (2005). 



